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Recently, we reported [M. Wagner et al. J. Mater. Res. 26, 1886 (2011)] transport measurements 
on the semiconducting intermetallic system Ruins and its substitution derivatives Ylulns-xAx {A 
= Sn, Zn). Higher values of the thermoelectric figure of merit [zT = 0.45) compared to the parent 
compound were achieved by chemical substitution. Here, using density functional theory based 
calculations, we report on the microscopic picture behind the measured phenomenon. We show in 
detail that the electronic structure of the substitution variants of the intermetallic system Rulns-^, A^, 
[A = Sn, Zn) changes in a rigid-band like fashion. This behavior makes possible the fine tuning of the 
substitution concentration to take advantage of the sharp peak-like features in the density of states 
of the semiconducting parent compound. Trends in the transport properties calculated using the 
semi-classical Boltzmann transport equations within the constant scattering time approximation are 
in good agreement with the former experimental results for RuIns-^^Sn^:. Based on the calculated 
thermopower for the p-doped systems, we reinvestigated the Zn-substituted derivative and obtained 
ZnO-free Kultis-xZiix. The new experimental results are consistent with the calculated trend in 
thermopower and yield large zT value of 0.8. 



I. INTRODUCTION 

Thermoelectric materials, capable of creating elec- 
tricity from waste heat are ideal in the search for sustain- 
able energy resources. The bulk of present day efficient 
thermoelectric materials have one common feature: all 
are heavily doped narrow band gap semiconductors with 
a carrier concentration of 10^^ to 10^^ cm~'^. In par- 
ticular, narrow band gap materials with complex crystal 
structure have been shown to possess excellent thermo- 
electric propertiesjii^ In general, the efficiency of a ther- 
moelectric material is measured by a dimensionless pa- 
rameter called the figure of merit, zT = S'^crT/ (nei+Hph), 
where S is the thermopower (Seebeck coefficient) of the 
material, T is temperature, a is the electrical conductiv- 
ity, Kei is the electronic part of the thermal conductivity 
and Kph is the lattice (phonon) contribution to the ther- 
mal conductivity. A zT value larger than 1 is impera- 
tive for successful applications. Since all the material- 
related parameters that determine zT are interlinked, 
tuning each parameter quasi-independently without ad- 
versely affecting the other is a tricky endeavor. 

Many of the binary compounds formed by transition- 
metal atoms with group III, IV and V elements have 
been shown to possess electronic band gaps. For ex- 
ample, FeSi, FeSb2 and FeGaa are all semiconducting 
with energy gaps of about 0.06 eV, 0.04 eV and 0.45 
eV respectively!^^— The hybridization of the narrow d 
bands of the transition metals with the rather broad 
p bands of the p-elements has been suggested to assist 
the formation of such band gapsi^ In particular, such 
hybridization also produces sharp features close to the 
Fermi level, which in turn has been shown to be quite 
beneficial for enhanced thermoelectric properties j^ii Con- 
sequently, large thermopower values of -1-500 /iV/K at 50 
K in FeSi,.^ -45 mV/K at 10 K in FeSbz,-^ and -350 ^lY /K 
at 300 K in FeGag^ have been evidenced. The compound 
Rulua is isotypic to FeGaa and was first reportectiS in 



1959 and later confirmed by HoUeck and co-worker&ii 
in 1964. Early resistivity measurements reported poor 
metallic conductivity in this system^ but later, density 
functional theory based calculations predicted semicon- 
ducting behavior Subsequent measurements have con- 
firmed the semiconducting nature of Ruins J^ii^ Studies 
on the thermoelectric properties of Rulna do not exist. 
Recently, we reported the results of the transport mea- 
surements of both Rulna and its substitution variants 
KuIn^-xAx {A = Sn, Zn).— The substitution variants 
exhibit relatively large Seebeck coefficients in a wide tem- 
perature range along with reduced thermal conductivity 
compared to the parent semiconductor Rulna , thereby in- 
corporating them in the family of potential thermoelec- 
tric materials. Transport measurements on the binary 
Rulua show semiconducting behavior with a band gap 
of ?a 0.45 eV for polycrystalline samplesJ^ The binary 
compound shows multi-band features exhibiting a large 
negative thermopower S = -363 //V/K at 308 K and a 
positive thermopower S = +262 /iV/K at 461 KJ^ This 
crossover behavior (electron transport at 308 K to hole 
transport at 461 K) is suppressed by chemical substitu- 
tion. The Sn-substituted samples show n-type behavior, 
while the Zn-substituted samples show p-type behavior, 
with a transition from semiconducting to metallic be- 
havior with increasing x. Chemical substitution reduces 
the thermal conductivity to 50% of the value of the bi- 
nary compound, resulting in a zT of 0.45 at 630 K for 
RuIn2.95Zno.05; an improvement by a factor of 7 over pure 
binary phase. The wide chemical versatility of this mate- 
rial, combined with the observance of good thermoelec- 
tric properties, suggests that further investigations are 
necessary to optimize and tune this material and identify 
similar versatile members (abundant and low-cost alter- 
natives) as viable candidates for future high-temperature 
thermoelectric applications. 

Here, firstly we present our results on the electronic 
and transport properties of Ruins and the substitution 
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FIG. 1. (Color online) Crystal structure of Ruins. Top panel 
shows the distorted cubes formed by In2 atoms, with Inl 
atoms sitting in the centers. Bottom panel shows the dis- 
torted trigonal prisms formed by In2 atoms with Ru atoms 
at its centers. Two neighboring trigonal prisms share a rect- 
angular face. The prisms and the cubes are interlinked along 
the c axis by the corner In2 atoms. 



variants using density functional theory (DFT) based cal- 
culations to provide a microscopic insight into the physics 
behind the measured phenomenon. The three theoreti- 
cal reports on the parent compound that exist in litera- 
ture have mainly focussed on evaluating the band struc- 
ture to obtain an estimate for the band gapJ^ii^iii In 
the present work, we focus on investigating the changes 
to the electronic structure as a function of substitution, 
as well as to probe the dependence of thermoelectric 
transport properties on the carrier concentration based 
on first-principles calculations. By doping a semicon- 
ductor, one creates a slightly different material and the 
changes induced by doping manifest themselves as a tran- 
formation of the band structure and, correspondingly, 
the electronic density of sates. Various physical mecha- 
nisms, including electron-electron interaction, impurity- 
band widening, band tailing, and screening effects be- 



tween the dopant and the host structure play a role in 
the altering of the band structurei^ Additionally, the 
dopant can either be randomly distributed or enter the 
host structure in an ordered fashion. These varied sce- 
narios can have a profound effect on the electronic struc- 
ture of the substituted variants, and the resulting band 
structure could either be a simple rigid-band shift of the 
semiconducting host compound to account for the addi- 
tional holes / electrons or a complex scenario with unequal 
band shifts of the intrinsic band edges. Tuning thermo- 
electric properties of the parent semiconductor by dop- 
ing, requires a well controlled change in the density of 
states of the doped variants, thereby accessing the sharp 
peak-like features close to the Fermi level of the host ma- 
terial. Consequently, it necessitates a careful study of 
the electronic structure of both the doped and un-doped 
compounds first and the subsequent calculation of the 
transport properties. In our work, we follow this two- 
step approach for Kvlii^-xAx {A = Sn, Zn). 

In our previous work, we observed small amounts of 
a ZnO impurity phase and non-reacted elemental Ru in 
the RuIn3_j;Zn2; samples The presence of these im- 
purity phases may result in the somewhat non-uniform 
behavior of the resistivity and thermopower values (as a 
function of x) for RuIna-ajZnj: in the temperature range 
300 K < T > 600 K. On the contrary, thermopower of 
the Sn-substituted samples reduced uniformly with in- 
creasing X. Based on our calculations of the transport 
properties of these materials (section V), we infer that 
the thermopower of both n- and p-doped systems should 
show a uniform behaviour. Consequently, we have suc- 
cessfully prepared a new batch of ZnO free RuIn3_^Zna; 
{x = 0.025, 0.050 and 0.100) applying starting materials 
of higher purity and measured their transport properties. 



II. CRYSTAL STRUCTURE 

Subsequent to the initial refinement of the crystal 
structure of Rulna within the space group PAn2 (CoGaa 
type), it has been shown that the this compound crys- 
tallizes in the higher symmetry of the tetragonal space 
group P4:2/mnm (no. 136), and belongs to the FeCaa 
type of structureJ^ The lattice parameters used in our 
calculations are a = 6.999 A and c = 7.246 A The unit 
cell of Ruln3 contains four formula units, with two crys- 
tallographically inequivalent In sites, Inl and In2 (Ru: 
[0.3451, 0.3451, 0]; Inl: [0, 0.5, 0]; In2: [0.1555, 0.1555, 
0.2622] )J^ The structure of Rulna can be visualized us- 
ing two basic building blocks in a complementary fashion: 
(i) each Inl atom sits in the center of a distorted cube 
formed by eight In2 atoms (Fig.[T] top panel); (ii) around 
each Ru atom six neighboring In2 atoms form a distorted 
trigonal prism, with two trigonal prisms sharing a com- 
mon face. Additionally, the trigonal prisms are inter- 
linked along the crystallographic z direction by sharing 
In2 corners (Fig.[T] bottom panel). Both building blocks 
share side faces and form a complex three-dimensional 
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FIG. 2. (Color online) Band structure (left panel) and 
total and site projected electronic density of states (right 
panel) of non-spin-polarized Ruins. The band structure 
is plotted along the standard high-symmetry directions (in 
units of (tt/o, 7r/b, TT, c)) of a tetragonal lattice: r(0,0,0) 
=^ X(0,l/2,0) ^ M(l/2,l/2,0) =^ Z(0,0,l/2) =^ 

R(0,l/2,l/2) ^ A(l/2,l/2,l/2). 



packing of polyhcdra, rather than a layered one. 



III. CALCULATIONAL DETAILS 

Non-spin-polarized DFT total energy and Kohn- 
Sham band-structure calculations were performed ap- 
plying the full-potential local-orbital code (version 
FPLO9.01.35), within the local density approximation 
(LDA)J^iS The Pcrdcw and Wang flavor^ of the ex- 
change correlation potential was chosen for the scalar 
relativistic calculations.— The calculations were carefully 
converged with respect to the number of k points, and 
finally a 24x24x24 regular /c-point mesh (1183 points 
in the irreducible Brillouin zone (BZ)) was used in the 
full BZ. The lattice parameters for undoped Rulna and 
as well as the doped variants^^ are taken from Ref.fTsl. 
It is well known that LDA tends to ovcrbind, resulting 
in lattice constants that are too small. Therefore, the 
calculations were carried out with experimental lattice 
constants, while the internal coordinates are relaxed for 
all the structures considered in this work. The substi- 
tutional derivatives RuIna-ajSn^: and RuIn3_xZnj. were 
modeled using three kinds of approaches: the simpler vir- 
tual crystal approximation (VGA), the fully disordered 
coherent potential approximation (CPA, calculated using 
FPLOS.OO)^, and the ordered supcrcell (SC) approach. 
Within VGA, a virtual atom with the desired number of 
electrons and protons is substituted at the indium wyck- 
off positions. We employed the VGA in the following 
way: the input of our full potential method are the nu- 
clear positions and charges and the basis functions (lo- 
cal orbitals). A virtual atom is constructed such that 



the number of effective valence electrons equals the one 
in the disordered alloy approximated by VGA. This is 
achieved by choosing an appropriate non-integer nuclear 
charge and a corresponding equal number of electrons 
for the disordered sites. The full potential scheme then 
solves the DFT equations for this unit cell with non- 
integer nuclear/electron charges, automatically adapting 
the basis functions and the potential to this nuclear con- 
figuration as it does for any other arrangement of "nat- 
ural" atoms. This choice of VGA reproduces the correct 
electron filling and often gives very satisfactory results 
compared to supercell and/or GPA calculationsi^i^ In 
GPA, the disorder of the dopant is taken into account 
explicitly and hence more reliable, but also more time 
consuming. Within GPA, the dopant atom retains its 
identity, but produces an effective medium that accounts 
for random disorder, resulting in the incoherence of the 
bands or more precisely, the spectral function A{k,Ltj). 
The incoherence arises because of the quantum mechan- 
ical averaging of the wavefunctions over the atomic po- 
tentials. In the SG approach, one of the indium posi- 
tion is occupied by the substituent Zn or Sn. Such a 
construction results in ordering of the substituent, but 
allows one to quantify specific effects that arise because 
of the periodicity that is imposed. Two supercells were 
constructed. The first was constructed by substitution of 
one In by Zn(Sn) within the tetragonal P4,2/mnm struc- 
ture unit cell of Ruins. Since one unit cell contains four 
formula units, this corresponds to 8.33% Zn(Sn) substi- 
tution and has no Zn-Zn(Sn-Sn) nearest-neighbor pairs. 
The second cell is a doubled cell along the c axis with 8 
formula units per unit cell. Again, one In was replaced 
by Zn(Sn), corresponding to 4.2% Zn(Sn) substitution. 
Like the smaller supercell, this cell also has no nearest- 
neighbor Zn-Zn(Sn-Sn) bonds. 

The transport properties were calculated using the 
scmiclassical Boltzmann transport theorj«2Sri2S within the 
constant scattering approximation as implemented in 
BoltzTraPi^S This approximation is based on the assump- 
tion that the scattering time r determining the electrical 
conductivity does not vary strongly with energy on the 
scale of kT. Additionally, no further assumptions are 
made on the dependence of r due to strong doping and 
temperature. This method has been successfully applied 
to many narrow band gap materials including clathrates 
and as well as to oxides 



IV. ELECTRONIC STRUCTURE 
A. Ruins 

GoUected in Fig. [2] are the LDA non-spin-polarized 
band structure and density of states (DOS) of the tetrag- 
onal Rulna. In the energy range displayed (-3.5 eV 
to 2 eV) there is a strong hybridization between the 
In 5p states and the Ru id states. We obtain a band 
gap value of 0.3 eV, consistent with previously published 
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FIG. 3. (Color online) Top panel: Calculated electronic DOS 
for the substitution derivative RuIn2.875Sno.125 and the binary 
compound Rulna as a reference. The vertical line at zero en- 
ergy denotes the Fermi level and has been set to the bottom 
of the conduction band of the parent compound for easy vi- 
sualization of the changes upon electron doping. The features 
close to the Fermi level are quite similar for VCA, CPA and 
the SC approach and change the DOS in a rigid-band like 
fashion compared to the parent compound. The difference 
between VCA and SC approach below -0.3 eV arise from the 
fact that, within SC approach the introduction of an actual 
Sn atom with a larger extent of the 5p orbitals increases the 
4d-5p band width. Disorder induced broadening is also ob- 
served for the CPA DOS. Bottom panel: The non-magnetic 
site projected band structure of RuIn2.875Sno.125 using the 
SC approach. For easy comparison with the band structure 
of the binary compound in Fig. (2] the SC approach bands 
have been unfolded into the smaller tetragonal unit cell. The 
bands crossing the Fermi level have both In and Sn character. 



results.— i^iii It is well known that the Kohn-Sham 
bands obtained from eonteniporary exehange eorrelation 
approximations underestimate the size of the semicon- 
ducting band gaps, hence the calculated band gap of 
0.3 cV is naturally slightly smaller than the experimen- 
tally observed value of 0.45 eV.— From the band struc- 
ture, we witness that Ruins has an indirect band gap 
with the top of the valence band occurring close to the 
A point along the R^=>A symmetry line and the bottom 
of the conduction band is close to the Z point along the 
r =^Tj line. It is worthwhile to note that the curvature 



of these two bands are similarly parabolic, resulting in 
similar values for the effective masses for the holes and 
electrons at low substitution concentrations. Our results 
concerning the bottom of the conduction band are consis- 
tent with the reports of Imai and co-workers^^ and with 
that of Mani and co-workersj^ On the contrary, calcu- 
lations by Bogdanov et alM using the TB-LMTO-ASA 
(tight binding-linearized muffin tin orbital-atomic sphere 
approximation) package find the rather flat band along 
M^=^> r as the minimum of the conduction band thereby 
yielding a larger effective mass of the holes as compared 
to electrons. Additionally, the bands are similarly dis- 
persive along the various high symmetry directions, con- 
sistent with the three dimensional nature of the build- 
ing blocks in the crystal structure. Our calculations are 
well in agreement with the recent experimental results on 
Ruing .i^ 



B. Ruing ajSria;: Electron doping 

The consequences of substituting Sn (group IV) for 
In predicted by the three calculational approaches are 
displayed in Fig.|3l The top panel compares the total 
DOS for the various approaches for the 4.2% substitu- 
tion scenario RuIn2.875Sno.125. Compared to In, Sn has 
one additional electron and acts as an electron donor in 
the substitution variant. For easy comparison, the DOS 
of the parent compound is also included in the figure. 
RuIn2.875Sno.125 is metallic compared to the parent com- 
pound which is semiconducting (see Fig. [2]). The shape of 
the DOS close to the Fermi edge is quite similar between 
VCA, CPA and the SC approach for RuIn2.875Sno.125. 
The slight dissimilarity in the DOS away from the Fermi 
edge arise from the underlying difFercnccs in the method- 
ologies. Explicit inclusion (as compared to VCA) of the 
Sn atom with the more extended 5p orbitals in the or- 
dered SC approach results in a larger band width of the 
Ru Ad - In 5p - Sn bp hybridized states. This is clearly 
discernible in the DOS below -0.3 eV. Similarly, disorder 
induced incoherence of the valence band is also noticed 
for the CPA derived DOS below -0.3 eV. As mentioned 
previously, there are two crystallographically inequiva- 
lent indium positions in this structure (Inl at 4c and In2 
at 8j). Hence, for each SC, two different calculations were 
performed by placing the substituent Sn once at the 4c 
position and once at the 8j position. The Sn atom prefer- 
entially occupies the 8j (In2) position, which is energeti- 
cally favorable by 20 meV per formula unit compared to 
the 4c (Inl) position. The bottom panel of Fig. |3] shows 
the non-magnetic band structure of RuIn2.875Sno.125 ob- 
tained using the SC approach. For a 4.2% substitution, 
the original unit cell of the binary compound has been 
doubled along the z axis, resulting in a smaller Brillouin 
zone (BZ). This has two consequences, firstly there are 
twice the number of bands and secondly all the bands are 
folded back into the small BZ of the SC. The resulting 
band picture gets quite complicated and hard to discern 
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meaningful information from it. To facilitate easy com- 
parison with the band structure of the binary compound, 
we have unfolded the bands according to the translational 
symmetry of the original unit cell using the approach de- 
scribed in Refs.ll^ and H^- The Ru and In character 
unfolded bands of the SC in Fig. [3] are comparable to 
the bands in Fig. [31 though most of the degeneracies arc 
lifted due to the presence of the Sn substitucnt. The size 
of the symbols in the band structure plot refers to the 
fc-resolved weight of all the orbitals of the chosen site. 
The larger the size of the symbol, the larger is the con- 
tribution of the site to the band structure. Based on 
this, we can conclude that the bands crossing the Fermi 
level have both In and Sn character, comprehensible due 
to the fact that Sn acts as an electron donor, making 
the system metallic and is only one atomic number away 
from In in the periodic table, and hence similar to In. 

Comparing the DOS of the binary Rulna with that of 
RuIn2.875Sno.125, we observe a rather rigid-band-likc shift 
in the electronic structure. Similar behavior is observed 
for other substitution concentrations (not shown here). 
This result is extremely interesting from a technological 
point of view. Mahan and Sofo^ showed that a narrow 
distribution of the energy of the electrons participating in 
the transport process is needed to maximize thermoelec- 
tric efficiency. The binary Ruins possesses sharp peak- 
like features close to the conduction band minimum and 
as well as near the valence band maximum. Electron dop- 
ing, as shown here using Sn as a substituent, results in a 
rigid-band shift of the DOS without adversely affecting 
the sharp peak-like features, thereby opening the possi- 
bility of fine-tuning the substitution to achieve maximum 
thermoelectric efficiency. This will be discussed in detail 
in section V. 



C. RuIna-j^Ziia;: Hole doping 

Following the same steps as described above for the 
Sn substitution, we have calculated the electronic struc- 
ture of the Zn substituted system as well. Collected 
in Fig.|3] are the DOS and unfolded band structure for 
the 4.2% substitution scenario RuIn2.875Zno.125. Sub- 
stitution of Zn introduces holes into the system making 
RuIn2.875Zno.125 metallic compared to the binary Ruins. 
Similar to the observations in the electron doped system, 
the DOS close to the Fermi edge in the hole doped variant 
is comparable between VCA, CPA and the SC approach. 
Moreover, the 3d orbitals of Zn are less extended and fully 
filled, resulting in a strong In character of the bands that 
cross the Fermi level. Nevertheless, the changes in the 
DOS close to the Fermi level are still rigid-band like, ben- 
eficial in tuning transport properties. Contrary to the Sn 
doping scenario, Zn atoms desirably occupy the 4c (Inl) 
position which is energetically favorable by 27 meV per 
formula unit compared to the 8j (In2) position. This 
difference in preferred site occupations between the two 
substituents can be construed as follows: In2 atoms form 
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FIG. 4. (Color online) Top panel: Calculated electronic DOS 
for the substitution derivative RuIn2.g75Zno.125 and the bi- 
nary Rulus for reference. The vertical line at zero energy 
denotes the Fermi level and has been set at the top of the 
valence band of the parent compound for easy visualization 
of the changes upon hole doping. VCA, CPA and the SC ap- 
proach result in a metallic ground state with similar features 
of the DOS close to the Fermi edge. Analogous to the elec- 
tron doping results, the changes in DOS for the hole doped 
system are again rigid-band like. In contrast to Sn doping 
(see caption of Fig. [3]), Zn 4d orbitals are not that extended 
and hence the valence band spectrum is quite similar for both 
VCA and SC approaches. Bottom panel: The non-magnetic 
site projected band structure of RuIn2.875Zno.125 using the 
SC approach. For easy comparison with the band structure 
of the parent compound in Fig.O the SC approach bands 
have been unfolded into the smaller tetragonal unit cell. The 
bands crossing the Fermi level have a strong In character. 



a distorted trigonal prism with a Ru atom at its center 
(refer Fig.[T]). Extended p orbitals of the anions are nec- 
essary to facilitate the hybridization with the narrow d 
bands of the transition metal cations at the center of the 
distorted trigonal prism. Hence, Sn atoms which possess 
extended and partially filled p orbitals preferentially oc- 
cupy the In2 position. Contrariwise, the transition metal 
Zn atom lacks extended p-like orbitals and hence prefer- 
entially occupies the Inl position, which is at the center 
of a distorted cube made of 8 neighboring In2 atoms. Our 
theoretical observation of preferential site occupation for 
the Zn atom is consistent with the experimental obser- 



6 



300 



200 



100 



i4 

I « 



-100 



-200 



-300 



1 ' 1 



AAAAAAAAAAAAAAAJ 
AAAAAAAAJ ' 



A** ^aA^I^aa** Ruin Zn 

*A*%^* ■ 0.010 A 

it--' ■ 0.025 A 



■ O.ObO A 

. ■ 0.100 A 

Ruin Sn o.l25 



300 



200 



400 

T(K) 



600 



800 



FIG. 5. (Color online) Thermopower S as a function of tem- 
perature for both Zn and Sn substitution variants Kulns^xAx. 
A negative sign in S represents n-type behavior and a positive 
sign in S represents p-type behavior. This figure is directly 
comparable to the experimental data from Ref.lTsl. 



vation in another isotypic system CoIn3_j.Zn3;j2^ Using 
X-ray diffraction and neutron powder diffraction exper- 
iments, the authors observed that the substitution of In 
by Zn takes place in an ordered fashion, producing col- 
ored variants of the parent compound Coins with the Zn 
atom entering exclusively the position corresponding to 
the cube centers. 

At this juncture, it should be noted that the procliv- 
ity of the substituent to certain indium sites is discerned 
only from the differences in the total energies of the vari- 
ous supercell configurations. No strongly distinguishable 
relevant features are noticeable in the DOS plots when 
Sn(Zn) occupy either Inl or In2 position and a similar 
rigid-band like shift of the DOS compared to the binary 
Rulna is noticed for all considered configurations. This 
observation is relevant for the experiments, where the fine 
tuning of the doped samples can be performed without 
the fear of adversely affecting the electronic structure of 
the final compound. Additionally, the metallic radii of 
In and Sn are similar («1.60A) and hence size effects 
do not play a significant role for the electron doped sys- 
tems. On the other hand, the metallic radius of Zn is 
about 13% smaller (Ril.40 A) than that of In and hence 
Zn is incorporated easily as a substitute for In in Ruins. 
Owing to the smaller metallic radius of the dopant Zn, 
there is a likelihood for local structural relaxation around 
the dopant which could have a detrimental effect on the 
thermoelectric transport by altering the electronic prop- 
erties of the host. We have investigated this in detail, by 
allowing for the local relaxation around the Zn atom in 
RuIn3_a;Zn2; with x = 0.125. In the un-relaxed structure, 
the Zn atom (occupying the Inl position) is surrounded 
by a distorted cube of In2 atoms (see FigH]) with four 
In2 neighbors at 3.16 A and another four In2 neighbors 
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FIG. 6. (Color online) The main panel shows the calculated 
thermopower S and the inset shows the electrical conductivity 
(a) relative to relaxation time (r) resolved along the crystal- 
lographic a (Sxx, {a/r)xx) and c (Szz, {a/T)zz) directions. 



at 3.26 A. The Ru-Zn bond length is 2.65 A. The struc- 
tural relaxation reduces the Zn-In2 bond length by just 
2%, while the Ru-Zn bond length on the other hand is 
slightly increased to 2.70 A. Nevertheless, these changes 
did not alter the electronic structure significantly, likely 
since the trigonal prisms centered around Ru and the 
cubes centered around Inl site are heavily distorted, al- 
ready in the undoped system. 



V. THERMOELECTRIC PROPERTIES 

A. Calculations 

Having an improved understanding of the electronic 
properties of Unlns-xAx {A = Sn, Zn) systems, we 
now calculate the transport coefficients using the semi- 
classical Boltzmann theory and the rigid-band approach. 
Such a procedure has been shown to successfully predict 
optimal doping levels in other thermoelectric materials, 
for example clathrates^^ and Sb2Te3i^ One main concern 
when calculating transport coefficients, is the underesti- 
mation of band gaps using the standard DFT functionals. 
Such an underestimation of band gaps manifests in the 
reduction of thermopower at higher temperatures due to 
bipolar conduction. To overcome this problem, we con- 
centrate on analyzing the trend in the calculated ther- 
mopower for various concentrations of electron and hole 
doping rather than quantifying them. This approach re- 
sults in a consistent scenario between experiment and 
theory as shown below. Firstly, as seen in Fig.[2J it 
is important to note that the valence and conduction 
band of the stoichiometric compound are parabolic. As 
shown above, both Sn and Zn substitution act in a rigid- 
band like fashion. Therefore, for small amounts of sub- 
stituents, the bands will continue to remain parabolic. 
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FIG. 7. (Color online) Calculated electrical conductivities (a) 
relative to relaxation time (r) as a function of temperature 
for both Zn and Sn substitution variants of Kuln^-xAx. 



FIG. 8. (Color online) Calculated power factor S^ct with re- 
spect to relaxation time r as a function of temperature for 
both Zn and Sn substitution variants of Knlns-xAx- 



Hence, experimentally the transport can be modeled us- 
ing simple parabolic band expressions to obtain the con- 
centration of charge carriers in the system. We also no- 
tice that close to the conduction band minimum, there 
are multiple heavy bands (i.e. bands with less disper- 
sion and hence larger effective mass), which is advan- 
tageous for thermoelectric transport. In contrast, the 
valence band maximum has a larger energy difference to 
the next available heavy band. Consequently, from just 
the view point of electronic structure, an n-type material 
would retain better thermoelectric transport coefficients 
than a p-type for low doping concentrations and low tem- 
peratures. 

The values of the calculated thermopower are not 
affected by the constant scattering time approximation 
used to calculate the various transport properties, since 
the expression for S is independent of t. This means that 
S is directly dependent on the electronic structure of the 
material. Collected in Fig. [5] is the temperature depen- 
dence of thermopower for various substitution concentra- 
tions of Zn and Sn in Ruins • Sn substitution introduces 
additional electrons into the system resulting in nega- 
tive values of thermopower. Correspondingly, we obtain 
positive values of thermopower for Zn substitution which 
introduces holes in the system. Due to the lack of Hall 
data on all the measured samples, the assignment of the 
substitution concentration for the plots here are done by 
comparing the experimental thermopower measurements 
with the calculated values. For all substitution concen- 
trations considered here, the magnitude of S increases 
steeply with temperature until 300 K. Above 300 K, in 
RuIn3_a;Sn2; jS*! slightly reduces or remains constant. On 
the other hand, for RuIna-ajZn^; IS"! above 300 K contin- 
ues to increase with a smaller slope. Below room tem- 
perature and for x < 0.050, \S\ for the electron doped 
system is 1.5 times that of the hole doped system, con- 



sistent with the electronic band structure. In contrast, 
with increased doping {x > 0.050) and temperatures (T 
> 300 K), hole doped RuIna-^jZua; shows increased ther- 
mopower compared to the electron doped RuIn3_2;Sn2;. 
The maximum values of 1 51 (larger than 200 /iV/K) are 
obtained for the smallest doping concentrations in both 
p- and n-type. Our results, including the absolute val- 
ues of S*, as well as the trend upon increased doping for 
the n-type, are consistent with the recent experimental 
measurements!^ A direct comparison of calculations for 
the p-type with the recently published experiments is in- 
hibited by the presence of impurity phase ZnO in the 
measured samples and hence no clear trend is discernible 
from the experiments Presently, stimulated by our cal- 
culational results, we have successfully synthesized a new 
batch of ZnO free Kuln^-xZiix samples and measured 
their thermoelectric properties (details described in sec- 
tion VI) . The trend in measured S is now consistent with 
our calculated results for the p-doped systems too (refer 
to Fig.IlOja)). 

Furthermore, transport measurements on the parent 
compound Rulna find a change in the sign of S from 
n-type to p-type around 350 K. Similar behavior was ob- 
served for the related compound RuGasi^^ The authors 
assigned the sign change in S to the shift of the valence 
band to an extrinsic region along with the presence of 
light holes in the system and proposed a two-band model, 
one each for electron and hole, respectively for evaluat- 
ing S. In this scenario, holes with higher mobilities can 
influence S and change the sign from negative to pos- 
itive valuesi^ However, modeling transport properties 
including the presence of extrinsic charge carriers within 
a bulk system is more involved and beyond the scope of 
the present investigation. 

As discussed earlier, notwithstanding the tetragonal 
symmetry of the crystal structure of Ruins, the polyhe- 
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FIG. 9. (Color online) Microstructure of RuIns-aiZni, samples: (a) x = 0.025, polarized light; (b) x = 0.050, polarized light, 
the surface is perpendicular to the pressure direction during the SPS treatment; (c) x — 0.050, polarized light, the surface is 
parallel to the pressure direction during the SPS treatment; (d) x = 0.050, bright light. 



dra of the basic building blocks form a three dimensional 
packing. Previously, measurements on single crystals of 
Ruln3 have found a somewhat weak anisotropic resistiv- 
ity along [110] and [001] orientations Collected in Fig.lH] 
are the values of S and cr/r along the crystallographic 
a and c directions. Consistent with the three dimen- 
sional packing of polyhedra, the anisotropy is moderate 
in the calculated transport coefficients, with the p-doped 
system displaying a larger anisotropy than that of the 
n-doped variant. This result can be qualitatively under- 
stood by analyzing the band structures shown in Figs. [3] 
and m For the n-doped system, the bands crossing the 
Fermi level disperse in a similar fashion along the various 
symmetry directions. On the contrary, for the p-doped 
system, the bands crossing the Fermi level are more dis- 
persive in the a — b plane (X=>M=> F) compared to 
the c diretion (F =^Z). 

Collected in Fig. [7] are the calculated electrical con- 
ductivities a relative to the relaxation time r for the same 
set of systems considered in Fig.[S] Note that a has a 
linear dependence to r. Assuming a constant r for the 
various substitution variants, the increasing conductiv- 
ity for either hole or electron doping reflects the fact that 
more charge carriers are present for the transport process. 
Due to the presence of a rather flat band with a larger 
effective mass close to the conduction band minimum, 
Sn substituted systems exhibit smaller electrical conduc- 
tivities compared to the Zn substituted systems at low 
temperatures and low substitution concentrations. The 



trend observed for the thermopower previously is thus 
different than that of the conductivity, which will result 
in an interesting trend for the power factor (S'^cr), the 
numerator in the expression for the thermoelectric figure 
of merit zT. To this end, we have plotted S^a relative to 
the relaxation time r in Fig.[8l For temperatures up to 
300 K, the power factors of the Sn substituted systems 
are enhanced compared to those of the Zn substituted 
systems. Beyond 300 K, the trend is reversed and the 
hole doped systems exhibit an enhanced power factor. A 
quantitative estimate for the dimensionless figure of merit 
zT from DFT, entails an accurate estimation of Kph, the 
lattice thermal conductivity. This involves solving the 
phonon Boltzmann transport equations exactly and re- 
quires an accurate description of the harmonic and anhar- 
monic interatomic forces, which is presently a formidable 
task and is left for a more extended study in the fu- 
ture. Additionally, disorder and nano-structuring effects 
that arise during experiments modify Kph in a not so 
straightforward way for a concise theoretical description. 
Nevertheless, clues provided by the transport coefficients 
that are directly dependent on the electronic structure 
have been used to successfully tune the thermoelectric 
properties of other known materials<2^^i^ In the case 
of Ruins _ J, 5*71 J, and RuIn3_j,Zna;, our calculations of the 
power factor establish a robust scenario for both p and n 
doped systems with improved values for certain doping 
levels. 
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FIG. 10. Measured thermoelectric properties of RuIns-jjZria; 
as a function of temperature, (a) Seebeck coefficient 5", (b) 
electrical resistivity p, (c) total thermal conductivity k, and 
(d) dimensional figure of merit zT. 



B. Experimental data for RuIn3_a;Zna; 

The chemical analysis of the commercially available 
starting materials (Ru, In and Zn) showed a measur- 
able amount of oxygen content only in the case of Ru: 
0.54 mass % of oxygen, which is presumably present in 
the form of Ru02. This oxide impurity is completely 
removed via hydrogen reduction at 900 °C for 5 hours 
in a home-made set up (O2 not detectable after H2 
reduction). After this hydrogen treatment, the sam- 
ples were prepared strictly in accord to the description 
in Ref.lTH. Microstructure analysis of RuIn2.975Zno.025 
(Fig.inti) show single-phase homogenous products with- 
out any ZnO impurity. Only in the case of RuIn2.95Zno.05 
and RuIn2.9QZno.10 small trace amounts of residual Ru 
were found (Fig.inji), which might be due to the defect oc- 
cupation of the Ru site. The grains, with respect to their 



X 


a [A) 


c (Aj 


0.025 


6.9988(1) 


7.2456(2) 


0.050 


6.9988(1) 


7.2396(1) 


0.100 


6.9977(2) 


7.2422(3) 



TABLE I. Refined lattice parameters for RuIns-iZn^ 



size, shape and orientation are quite similar for all the 
samples discussed here. Additionally, the microstructure 
is independent from the substitution concentration of Zn, 
and the direction of pressure during the spark plasma 
sintering (SPS) experiments (Fig.|9)D and c). The X-ray 
diffraction patterns of the RuIn3_a;Zna: samples were in- 
dexed with the PA2/mnm (no. 136) space group. Com- 
pared to the previously published results^ the a lattice 
parameter remains unchanged, while the c lattice param- 
eter is slightly reduced (Table. I) . For composition of the 
samples, see Ref.[3^ Because of the absence of ZnO im- 
purities, all of the Zn has been incorporated into the 
lattice, resulting in a pronounced reduction of the unit 
cell in the c direction. 

The transport-related quantities and the resultant 
figure of merit zT are collected in Fig.[TUl The tem- 
perature dependence of the electrical resistivity p shows 
metallic behavior for all the substitution derivatives, with 
increasing resistivity for increasing temperature. The ab- 
solute values of p shown here are reduced in comparison 
to our previously published resultsi^ due to the absence 
of ZnO impurity phase, which is an insulator with a large 
band gap of more than 3 eV. The measured thermopower 
S is positive (p-type) and is reducing with increasing x. 
At 350 K, the absolute values of Sexp change from 162 to 
150 /iV/K for 0.025 < a; > 0.100, compared to our calcu- 
lations where Stheo change from 170 to 115 pV/K. For all 
the samples considered here, S increases with increasing 
temperature and reaches a maximum around 600 K. This 
trend in 5 as a function of x and temperature are now 
consistent with our above mentioned theoretical predic- 
tions (Fig. [5]) . The thermal conductivity k is only slightly 
increased in comparison to Ref.fisl Previously, the pres- 
ence of ZnO impurities could have acted as additional 
scatterers to reduce k, which is presently removed from 
our samples. The temperature dependence of the dimen- 
sionless figure of merit zT increases strongly, reaching a 
value of 0.8 at 620 K for Ruins ,975Zno.o25- Albeit the 
trend in S is congruent between experiment and theory, 
the same cannot be inferred for the other transport re- 
lated quantities and hence zT, because of the complex 
depedcncy relations between these parameters. While 
the expression for S is free of parameters, p, k and hence 
zT are dependent on the scattering time r. Our calcu- 
lations which are based on the constant scattering time 
approximation cannot be expected to yield more realistic 
results. Nonetheless, we have demonstrated that insights 
from band structure calculations can be quite helpful in 
fine tuning thermoelectric properties. 
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VI. CONCLUSIONS 



We have analyzed in detail the electronic structure 
of the hole- and electron-doped Ruins intermetallic com- 
pound. Wc have used various methods to model the ef- 
fects of doping, starting from the simple VGA descrip- 
tion to the time consuming and the more reliable CPA 
approach. Effects of ordering of the dopant was also 
considered using the SC approach. Comparing the re- 
sults of VCA, CPA and the SC approach, we can con- 
clusively determine that the dopants introduce charge 
carriers and change the electronic structure in a rigid- 
band-like fashion. This scenario allows for the advanta- 
geous use of the sharp peak like features in the DOS to 
obtain good thermoelectric properties. Transport coeffi- 
cients (thermopower and power factor) calculated using 
the semi-classical Boltzmann transport equations, are in 
good agreement with the experimental results. Based 
on our calculations, we can conclude that low substitu- 
tion levels are advantageous in obtaining large values of 
thermopower. Both n- and p-doped systems have similar 
transport properties. In particular, based on our calcu- 
lations, for T > 400 K, p doped samples with x > 0.125 



show considerable potential for improving the figure of 
merit zT. Similarly, for T < 500 K, n doped systems 
with an optimal substitution in the range 0.050 < a; > 
0.100 also show great thermoelectric potential. A zT sa 
0.8 has been obtained for RuIn2.975Zno.025, while the first 
trial of experiments were able to only reach a zT « 0.1 in 
the n-doped samplesJ^ Additional experimental investi- 
gation exploring higher p-dopant concentrations as well 
as fine tuning the n-doped systems are desirable to probe 
the possibility of further enhancing zT. Kulns-x-^x = 
Sn, Zn) systems with their robust electronic properties, 
containing low-toxic combination of elements provide - 
although expensive - an ideal play ground for experimen- 
talists to investigate in detail and tune the thermoelec- 
tric properties as a function of doping. Our combined 
results of theoretical calculations and experimental work 
has clearly demonstrated that narrow-band-gap semicon- 
ductors formed by transition-metal atoms with group III, 
IV and V elements have a large potential in thermoelec- 
tric applications. The presented study will, hopefully, 
also stimulate the search for new better and cheaper sys- 
tems in this family of compounds whose properties follow 
the same underlying physics as Kuln^^x-'^x = Sn, Zn). 
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